Tails of the Density of States in a Random Magnetic Field 
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We study the tails of the density of states of fermions subject to a random magnetic field with 
non-zero mean with the Optimum Fluctuation Method (OFM). Closer to the centres of the Landau 
levels, the density of states is found to be Gaussian, whereas the energy dependence is non-analytic 
near the lower bound of the spectrum. 
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The problem of a charged quantum particle con- 
strained to move in a two dimensional (2D) static random 
magnetic field (RMF) has attracted considerable theo- 
retical and experimental interest in the past few years. 
The model plays an important role within the compos- 
ite fermion picture of the fractional quantum Hall ef- 
fect 0. Furthermore, it is supposed to describe states 
with spin-charge separation in high-Tc superconductors 
• It is also relevant to the understanding of the proper- 
ties of a two dimensional electron gas (2DEG) in lattice- 
mismatched InAs/InGaAs heterostructures in magnetic 
fields 0. In the latter systems the electron gas is non- 
planar due to the lattice-mismatched epitaxial growth. 
When a uniform magnetic field B is applied, the elec- 
trons experience an effective inhomogeneous field per- 
pendicular to the non- planar 2DEG [2j. In addition, a 
static RMF in 2D inversion layers can be experimentally 
realized in several ways. One possibility is to use a type- 
II superconductor with a disordered Abrikosov flux lat- 
tice in an external magnetic field as the substrate for the 
2DEG 3 . Alternatively, a magnetically active substrate 
such as a demagnetized ferromagnet with randomly ori- 
ented magnetic domains may be used :5] . Recently, static 
RMFs in 2DEGs were created by applying strong mag- 
netic fields parallel to GaAs Hall-bars decorated with ran- 
domly patterned magnetic films 

The most fundamental quantity for understanding the 
electonic properties of a random system is the density of 
energy levels. The standard method to estimate the den- 
sity of states (DOS) is to calculate the imaginary part 
of the trace of the single-particle Green function by di- 
agrammatic techniques. However, this approach fails in 
the tails of an energy band where multiple scattering up 
to infinite order has to be considered in order to take 
into account correctly the effect of localisation of elec- 
trons. Also, numerical approaches are bound to fail in 
the asymptotic tails since here the eigenstates are deter- 
mined by rare statistical fluctuations of the randomness. 
Moreover, in the case of RMF, the perturbative approach 
is also fundamentally problematic since one has to deal 
with the non-diagonal part of the Green function, which 
is not gauge invariant. In addition, the calculation of 
the__Green function is plagued by infrared divergencies 
ITol Hil l that are due to the long-range nature of 



the correlations of the vector potential, even if the spa- 
tial correlations in the RMF are short-ranged. It has 
been suggested that these divergencies are due to the 
non-gauge-invariance of the Green function and there- 
fore unphysical 0, although^;ecently, a physical inter- 
pretation has been proposed [O . In order to avoid such 
difficulties, E. Altshuler et al. calculated the DOS 

of a charge in a RMF using the semiclassical approxi- 
mation. This is valid when the energy E is much larger 
than the cyclotron energy hujc corresponding to the mean 
magnetic field B. Also, it should exceed F, the disorder 
induced width of the Landau Levels (LL). A field theo- 
retical approach has been used to determine the DOS in 
a RMF with zero mean value near the band edge Q . The 
tail of the DOS in a system of randomly distributed flux 
tubes of fixed strength was considered 0. There are 
also several numerical studies of th e sp ectrum with dif- 
ferent mean and correlation lengths [ij] . Recently, math- 
ematically rigorous results have been obtained ^15, d61] . 
In particular |0, upper and lower bounds for the loga- 
rithm of the integrated DOS near _B = of some simple 
Gaussian RMFs with zero mean values have been esti- 
mated. For RMFs with non-zero mean value, the limit 
when E is smaller than hujc and, more generally, the tails 
of the lower Landau bands have not been considered an- 
alytically so far. 

It is our purpose to provide nonperturbative results for 
the DOS in the tails of the lower Landau bands, as broad- 
ened by a static RMF. We show that the Optimum Fluc- 
tuation Method (OFM) 0,E3, being non-perturbative 
and free from divergencies, can be extended to treat this 
RMF problem. We consider non-interacting fermions in 
a RMF with non-zero mean, B + b{r), with b{r) Gaussian 
distributed. The Hamiltonian 
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has a sharp lower bound of the energy spectrum at E — 0. 
We concentrate on the energy region near the first Lan- 
dau band. Since the OFM is especially designed to grasp 
rare fluctuations, it allows us to calculate the energy- 
and i?-dependence of the leading terms of the DOS. The 
correlation function of the RMF is assumed as 



(6(r)fe(r'))-/3(|r-r'|). 



(2) 
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with /3(|r — r'l) ^ as |r — r'| oo. The random 
field is ergodic, i.e. the correlations between different re- 
gions decay to zero with increasing distance. We also 
assume that (3{r) is characterized by a single scale Tc, 
the correlation length of the RMF. The probability den- 
sity for a specific realization 6(r) with correlator |5J) is 
P[6(r)] = A/'exp{— S'[6(r)]}, where JV is the normaliza- 
tion constant and S'[6(r)], the action of the RMF, is 

S[b{r)] ^ J b[v)p-^[v - r')&(r')drdr'. (3) 

The kernel /3~^(r — r') has the property 

j dr'p-\r - r')/3(r' - r") = <5(r - r"). (4) 

For a (5-correlated RMF, the action reduces to 5[6(r)] = 
l//3o J b^{r)dr. The configurationally averaged DOS is 

piE) = J Db{v)P[b{v)]p{E- [6(r)]). (5) 

In the low energy tail of the first Landau band states 
are expected to be localized near strong, exponentially 
rare fluctuations of 6(r). Therefore, the average in Eq. 
(O over all configurations, yielding states at energy E is 
dominated by the most probable realization of 6(r). The 
functional integral in Eq. ||SJ| can be evaluated using the 
saddle point approximation. Furthermore, only fluctua- 
tions in which the lowest level Eq is equal to E have to 
be taken into account since a configuration in which E 
corresponds to an excited level is less probable. Within 
logarithmic accuracy [Ulil El El El El , 

-\np{E) ^ ^\nS[b{r)]\E,,[b(r)]=E- (6) 

f)(r) 

Before we present a rigorous treatment, let us first 
give an intuitive estimate for the DOS closer to the 
band center, AE <^ huc, where E = fiiOc — AE. For 
short-range RMF, rc < Ib, with Ig = {hc/eBfl'^ the 
magnetic length related to B, an optimal configuration 
near the band centre is likely to be a circular magnetic 
well with depth A6 ^ B. The corresponding action is 
S « 7rA6^i?^//3o, where R is the radius of the well. The 
ground state energy Eq of a charged particle in a circu- 
lar magnetic well of radius where the magnetic field 
is different from the constant magnetic field B, is known 
[23| . The radii i?opt of the optimal wells with lowest ac- 
tion S'min(A£'), are proportional to Ib and independent 
on Ai?. Moreover, the depths A6opt <^ B oi these wells 
are proportional to AE and do not depend on B, 

Ropt ^ 1.6 Ib cx AE cx Afeopt- (7) 

The action of the optimum fiuctuation with Eq — E is 
then 



From ® and ((SJ we obtain 

p{E) ~ exp (-AE'/TIq) , (9) 

with r^^o = CihefiQ^'^ / {mcclB)^ where a ^ 1.5 • 10~^ is a 
numerical factor. The variance of the Gaussian is thus 
proportional to B. Our simple arguments are expected 
to be valid when the energetic distance from the center 
of the lowest LL fulfils F^^o ^ <^ hwc- A completely 
analogous argument holds for the right tail of the first LL, 
near the band centre. In this case the optimal fiuctua- 
tions are magnetic circular humps with height A& <C B 
and radius R (x Ib and the leading exponential term of 
the DOS in the right tail shows the same dependence on 
the energy shift and B a.s 

With long-range RMFs, the analysis of the DOS near 
the band centre is simpler: the localization radius of 
a typical state of the order of Ib is much shorter than 
the radius of an optimal potential well. The correlation 
length Tc, and the energy E of such a state is, in lead- 
ing order, equal to the first LL energy in the total field 
B — Ab. The energy shift is thus proportional to Ab, and 
the RMF acts exactly like a random electrostatic poten- 
tial. Since for long-range RMFs the radius of the well is 
the largest length scale, the probability distribution Q 
can be approximated as 

P[Ab] - exp (-A5V/3(0)) . (10) 

Hence, 

AE'^ 

_lnp(iJ)^ (11) 
^ 

with Fq = fte/3(0)^/^/TOeC, and the exponent of the DOS 
does not depend on B, as long as the inequality Ib ^ 
is fulfilled. Equation l(TT|l is valid if Fq < AE ^ hujc- 
Similar considerations are expected to yield a Gaussian 
DOS also in the tails of higher Landau bands, in the 
regions F„ ^ — (n + l/2)hajc\ ^ huJc, where F„ is the 
width of the n-th LL. In these regions, the DOS resembles 
the one of independent charged particles in a Gaussian 
electrostatic potential |2l l2l |. 

Due to the sharp band edge at i? = the DOS 
is expected to approach zero more rapidly for energies 
E <^ hujc- For short-range RMF, states with arbitrarily 
small energies can be obtained when they are localized in 
regions of area A, inside which & « and outside which 
h K, B. The action of these fiuctuations is 

S - AB'^/pQ, (12) 

and the ground state energy scales like the one in a po- 
tential weU, E ~ r?l2m^A. Thus, XjAo^E and the DOS 
becomes a non-analytic function of E, 

p{E) ^ exp {-KqB'^IPqE) , (13) 



(8) 



with Kq ~ 'nti} jlnii,. The above picture is analogous 
to the argument used by Lifshitz to estimate the tail of 
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the DOS of a particle subject to a Poissonian random 
potential generated by short-r ang e, repulsive impurities 
in zero magnetic field [TM l2(l l2,'il | . The argument holds 
for a Poissonian distribution of magnetic fluxes, too [T3l| . 

Intuitively, for large correlation length, Vc > Ib, there 
is a right neighbourhood of = such that each corre- 
sponding optimum well has a radius much larger than 
the correlation length, R Vc- Equivalently, E <C 

/{mf,r1). The larger the correlation length, the closer 
to the band edge the energies of the optimal states must 
be in order to fulfil this inequality. Defining 



(14) 



the action of 6(r) is still given by Eq. H12|l and the DOS 
by Eq. (fT!^ . Hence, for £' — > 0, the DOS becomes inde- 
pendent of the correlation length r^. 

In order to obtain exact expressions for the DOS, we 
now derive the variational equations which determine the 
shape of the optimal fluctuations and wave functions [Tol 
l2ll |. According to ||HJ), we must search for the maximum of 
the probability distribution Eq. Q under the constraint 
Eq = E. For the weaker constraint 



det{iJ[fe(r)] -E}^{), 



(15) 



or, equivalently. En = E for some energy level En, the 
optimum fluctuation 6(r) of the RMF must satisfy 



l3-^{s-s')h{s')ds' 



Sdet{H[b{r)] - E} 



Sb{s) 



b=b 



0, 

(16) 



where /i is a Lagrange multiplier. Using 

det{H[b{r) - E} = exp(tr ln{i? [6(r)] - E}), (17) 

and assuming that the ground state energy £'o[&(r)] is 
equal to E, we find 



p-\s~s')b{s')ds' + ^'{E) 



SEoHr)] 



5b{s) 



(18) 



b=b 



with ^'{E) = A*ni^i(^n ^ the Coulomb gauge, 

we can write the Hamiltonian as a function of 6(r) and 
calculate SEQ/6b{s). The variational equation (|18|l yields 



b{s) = -fi'{E)^ I cfs'pis-s') I d^r jo-a<i>(r-s') (19) 



where 



a$(r — s) = 



1 z X (r — s) 
2^ |r-sP 



(20) 



and jo = ^'on^'o/2mo is the ground state current, 11 is 
the kinetic momentum of the particle. Eq. ((T^ , together 
with 
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\ 


~ 1/E 











0.0 



0.1 



0.2 



0.3 



0.4 



0.5 



E 



FIG. 1: Behavior of the effective action S{E) = — In p{E) in 
a short-range RMF, Vc — 0.11b- Energies are in units of hcuc 
and Fo ~ 10~'^fkJc, S is in arbitrary units. 



determines the optimal magnetic field 6(r) and ground 
state wave function 'l'o(r)- 

We have solved equations H19|) and (|21|l iteratively both 
in the case ^ Ib and in the case Vc ^ Ib for RMFs 
with Gaussian correlators 



P{r - r') = 



27rr; 



■ exp 



2r2 



(22) 



Since the RMF distribution function is rotationally in- 
variant, circular symmetry is assumed. 

For short-range RMFs the result for the DOS p{E) is 
shown in Fig. ^ as function of energy E, where S{E) = 
— \np{E). At the band edge we find indeed that the ac- 
tion of the optimum fluctuation shows the characteristic 
behavior as a function of the energy E, S ^ E^^ while 
closer to the band center it changes to S* oc AE^. This 
reflects the physical origins of the corresponding typical 
wave functions, and is consistent with the qualitative ar- 
guments, given above. 

Fig. 13 shows the optimal fluctuations near the band 
centre and the band edge for short-range fluctuations. 
Near the centre, the typical fluctuations are shallow wells, 
compared to B, with relatively steep walls. In the case of 
long-range flelds, near the band centre, the radius of the 
ground state is much smaller than the size of the well. 

In experiments some kind of random electrostatic po- 
tential is always present. Let us assume that the RMF 
and the random potential (RP) are independent random 
quantities and that they are both long-ranged. For a 
weak RP, W{Oy/^ < heP{Oy^^ /(mc) (where W{r) ^ 
{V{0)V{r))), the RMF is dominant in the tail at posi- 
tive energies except for a narrow region close to E = 0. 
At E hujc the action of an optimal well of the RP 
is Srp ^ {E- hjjj2f/2W{Q) - h^ujl/m{Q), whereas 
Srmf ~ n^Lol/{2mErl6ul), where Suu^ = ePiOf^ /mc 
and Tc is the correlation length of the RMF. The RMF 
will dominate on the RP if Srmf < Srp; therefore. 
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FIG. 2: The optimum magnetic wells b{r) (solid line) and the 
corresponding ground state wave functions ^'(r) (dashed line) 
for correlation length rc — 0.11b- Left: closer to the centre 
of the band, energy shift AE = 3.6 ■ 10~^huJc- Right: near 
the band edge; ground state energy _E = 1.8 • r is in 

units of Ib and 6(r) in units of B. 



the exponent of the DOS will be proportional to 1/E if 
1 < i?V^? < h^Suj^/W{0), with i?2 ^ h^/mE. Hence, 
the results presented in this paper break down at en- 
ergy Ec ~ W{0)/mr'^Suj'^. For large, negative energies, 
E — > —oo, the RMF becomes irrelevant and the DOS is 
purely classical, -lnp(£') = E^/2W{0). 

In conclusion, we have determined the density of states 
of a charged particle in a spatially correlated, randomly 
varying magnetic field with non-vanishing average. The 
latter provides Landau levels which are broadened into 
Landau bands by the RMF, well separated for sufficiently 
small disorder. In the regions of the tails of the Landau 



bands which are accessible neither by perturbative multi- 
ple scattering expansions nor by numerical calculations, 
we have found that the average DOS is determined by 
the typical configuration of the magnetic field. This is 
reflected in the energy dependence of the effective ac- 
tion, S{E), and the fact that the latter is proportional 
to the logarithm of the DOS, p{E), which is found to 
be asymptotically singular at the lower bound of the en- 
ergy spectrum and becomes quadratic as a function of 
the energy closer towards the centre of the band. In or- 
der to determine the pre-exponential factor of the DOS, 
one integrates over the fluctuations of the magnetic field 
around the saddle point configuration satisfying Eq. l(T^ 
[21I l2fil | . As a final remark, we want to discuss briefly 
the relevance of our work to the CF description of the 
Fractional Quantum Hall Effect. Within this model, elec- 
trons are replaced by fermions experiencing a fictitious 
magnetic field proportional to the particle density in ad- 
dition to the external one. In the presence of a random 
impurity potential, at the mean field level, the particle 
density is spatially inhomogeneous due to screening and 
the fictitious magnetic field has thus a spatially stochastic 
component 1]. However, since the random impurity po- 
tential and the RMF are not independent random quan- 
tities, our theory cannot be straightforwardly applied to 
the CF model. This issue is subject to future work and 
will be published elsewhere. 
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